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Solving null-geodesic equations, behavior of angular diameter distances is studied in 
inhomogeneous cosmological models, which are given by performing JV-body simulations 
^ . with the CDM spectrum. The distances depend on the separation angle 6 of ray pairs, 

the mass m and the radius r s of particles consisting of galaxies and dark matter balls, 
' and cosmological model parameters. The calculated distances are compared with the Dyer- 

Roeder angular diameter distance with the dumpiness parameter a (=0 — 1), and for each 
7-H ■ ray pair the corresponding a is determined. Shooting many ray pairs, we derive statistical 

^ ■ quantities such as the average value (a), the dispersion (a a ), and the distribution of a. 

|>. ; It is found in four cosmological models with {O , A ) = (1, 0), (0.2, 0), (0.4, 0), (0.2, 0.8) 

for the particle parameters m ~ 2 x 10 11 Mq and r B — (10 — 40)/i _1 kpc that (1) a is nearly 
equal to 1 or the Friedmann distance is best fitted, (2) a a decreases with the increase of 
\Q ' the redshift z and radius r s , (3) for 9 >> 1 arcsec, <j a is very small, but for 8 < 1 arcsec it 

is so large that the use of only the the Friedmann distance (a = 1) may cause some errors 
OO . for quantitative analysis of cosmological lensing, (4) the distribution of a is symmetric and 

■ asymmetric for a a < 0.5 and > 0.8, respectively, and (5) a a is smallest and largest in models 

with (1, 0) and (0.2, 0), respectively. The cosmological constant has the role of decreasing 

'—^ ' §1. Introduction 

c/3 ■ 

Recently many kinds of lens phenomena are observed in strong and weak forms of 
deformed images of quasars and galaxies due to lenses of galaxies, clusters of galaxies, 
^ . and large-scale matter distribution. Lensing is being used as a new observational 

\ tool for remote dark objects. 

For the analysis of such lens phenomena, angular diameter distances are often 
used because of their convenience, but their definition is not unique. The most 
representative ones are the Friedmann distance and Dyer-Roeder distance, which 
are given as the distances in the homogeneous and smooth matter distribution and 
the distance in low-density regions in the inhomogeneous matter distribution, respec- 
tively. 2 ) The behavior of these distances has been discussed in various papers. 3 ) " 6 ) 
It is important to clarify which angular diameter distance is most appropriate in the 
observation of our Universe. This problem depends on the separation angles of ray 
pairs, the mass m and radius r s of lens objects, like galaxies and dark matter balls, 
and cosmological model parameters (/2o,Ao). 

In this paper we first use for this purpose the inhomogeneous models which were 
produced in previous papers, 7 )' 8 ) using TV-body simulations, and we then consider 
many ray pairs which are made by solving null-geodesic equations in the above 
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mentioned models and from which the average properties of the angular diameter 
distances are derived, taking into account the values of 0,m,r s ,f2o and Ao- In §2 
the behavior of the Dyer-Roeder angular diameter distance with the dumpiness 
parameter a in various background models is given explicitly. In §3 the angular 
diameter distances in inhomogeneous models produced by the numerical simulations 
are derived. By comparing those with the Dyer-Roeder angular diameter distance, 
the value of a for each ray pair is determined, and its statistical behavior is presented. 
In §4 concluding remarks are given. 

§2. Dyer-Roeder angular diameter distance 
with the dumpiness parameter a 

The Dyer-Roeder angular diameter distance D A (z) satisfies 

^D A + ^(l + z) 5 af2 D A = 0, (2-1) 

where a is the dumpiness parameter and v is the affine parameter. 3 ) " 6 ) v is related 
to z by 

dz/dv = (1 + z) 2 [(f2 z + 1)(1 + zf - \qz{2 + z)} 1 ' 2 . (2-2) 

The cases a = and 1 represent the limiting distances in the empty region and the 
Friedmann distance in the homogeneous region. 

As we consider the rays received by an observer at the present epoch, we have 
the conditions 

D A (0) = 0, (dD A /dz) z=0 = c/H , (2-3) 

where Hq (= 100/i km s _1 Mpc -1 ) is the Hubble constant. 

Here we consider four models with (^ ,A ) = (1, 0), (0.1, 0), (0.2, 0) and (0.2,0.8). 
The above equations are easily solved in each model and we obtain the z-dependence 
of D A (z) for a = and 1 and the a-dependence of D A {z) for z = 0.5, 1, • • • , 5, which 
are shown in Figs. 1, 2 and Figs. 3 ~ 6, respectively. These dependences are used 
to interpret the results in the next section. 

§3. Angular diameter distances in inhomogeneous models produced by 

A^-body simulations 

For the quantitative analysis of image deformation in the lensing phenomena, 
we treated in previous papers various clumpy Friedmann models (cf. Refs. 7) and 
8), and also Refs. 9) and 10)). In the present paper we use the same models for 
the statistical analysis of the angular diameter distances. In them we considered 
inhomogeneities in periodic boxes including N = 32 3 particles, whose distribution 
was produced by performing the iV-body simulation in the background models S, 
01, 02, and L with (f2 ,X ) = (1, 0), (0.2, 0), (0.4, 0), and (0.2,0.8), respectively. For 
the simulation we used Suto's tree code 11 ) during the time interval < z < 5. The 
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Fig. 1. The z dependence of the angular di- Fig. 2. The z dependence of the angular diam- 
ameter distance Da for a = 1. Models eter distance Da for a — 0. Lines have the 

S(1,0), 01(0.2,0), 02(0.4,0) and L(0.2,0.8) same meaning as in Fig. 1. 

are denoted by solid, dotted, short dash, 
long dash lines, respectively. 

initial conditions were given in the CDM spectrum with the power n = 1 and the 
dispersion <rg = 0.94, using Bertschinger's software COSMICS, 12 ) where the Hubble 
parameter was specified as h = 0.5,0.7,0.6 and 0.7 for model S, 01, 02 and L, 
respectively. 

The present sizes of the periodic boxes and the particle mass are 

L = (32.5, 50, 39.7, 50)/i~ 1 Mpc, (3-1) 

and 

m (= pboL 3 /N) = (2.90, 2.11, 2.11, 2.11) x HF/i^Mq (3-2) 

for models S, 01, 02 and L, respectively, where pbo = 3Hq 2 [2o/(8ttG). 

In all models, particles consist of compact lens objects which are galaxies and 
dark matter balls with the same mass m and the same constant physical radius r s . 
The representative value of r s is 20/i _1 kpc, but for comparison the cases of r s = 
10/j -1 kpc and 40/i -1 kpc are also considered. In a previous paper 7 ) it was assumed 
that in model S particles consist of compact lens objects (20%) with r s = 20h~ 1 
kpc and clouds (80%) with the same mass m and r s > 20h~ 1 kpc. In this paper we 
take the model S including only compact lens objects with various r s . Model S with 
larger r s corresponds effectively to the case including clouds with larger r s . 
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Fig. 3. The a dependence of the angular di- 
ameter distance Da in the model S with 
(i?o,Ao) = (1,0). The epochs with z = 
0.5,1,2,3,4, and 5 are denoted by solid, 
dotted, short dashed, long dashed lines, re- 
spectively. 



Fig. 4. The a dependence of the angular di- 
ameter distance Da in the model Ol with 
(0.2, 0). Lines have the same meaning as 
in Fig. 3. 



The line-element of inhomogeneous models is expressed as 



ds' 



1 



= -(1 + 2(p/c z )c z dt z + (1 - 2(p/c z )a z (t){dx) z /[l + -K{x) 



2i2 



(3-3) 



where ip is the gravitational potential connected with the matter distribution through 
the Poisson equation, and K is the signature (± or 0) of the spatial curvature. The 
physical radius r s is expressed using the comoving radius x s as r s = a{t)x s and it is 
taken into consideration in the form of a softening parameter (r s ) in the gravitational 
force as l/[a(t)(x — x n ) 2 + (r s ) 2 ], where x n is the position vector of n-th particle. 

Let us consider a pair of rays received by the observer with the separation angle 9. 
By solving null-geodesic equations, which were obtained in the previous papers, 7 ^' 8 - ) 
the interval of the two rays at any epoch can be derived. If (Ax)± is the component 
of the deviation vector perpendicular to the central direction of the rays, the angular 
diameter distance Da is defined as 



D A = a(t)(Ax) 



_[1 + \k{x 



where the factor (1 — 2(p/c 2 ) is neglected, because \y>/c 2 
expression can be rewritten by use of y % (= a$x l /Rq) as 

Rq 



2 rv#, (3-4) 

<< 1 locally. The above 



D A = 



(l + z)F 



(3-5) 
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Fig. 5. The a dependence of the angular di- 
ameter distance Da in the model 02 with 
(0.4, 0). Lines have the same meaning as 
in Fig. 3. 



Fig. 6. The a dependence of the angular di- 
ameter distance Da in the model L with 
(0.2, 0.8). Lines have the same meaning as 
in Fig. 3. 



where F = 1 — \(R H /c) 2 (l - f2 - X )(yf 



a 



o(t ), and R = Lq/N 1 / 3 . 



Behavior of the angular diameter distance depends on the separation angle 9. In 
the cosmological observation for angular sizes of compact radio sources by Kellerman, 
13 )> 14 ) it was found to be on the order of milli-arcseconds, or 9 ~ 0.005 arcsec. 
Accordingly we adopt 9 = 0.005 arcsec as a representative value and for comparison 
9 = 0.1 and 20 arcsec as other representative values. 

In the present lensing simulation we considered five parameter sets (r s ,9), A. 
(20/t _1 kpc, 0.005 arcsec), B. (40/t _1 kpc, 0.005 arcsec), C. (10/i _1 kpc, 0.005 arcsec), 
D. (20/t _1 kpc, 0.1 arcsec) and E. (20/i _1 kpc, 20 arcsec), and performed the ray- 
shooting of 500 ray pairs for each parameter set in the four models. At the six 
epochs z = 0.5, 1, 2, 3, 4 and 5, we compared the calculated distances with the Dyer- 
Roeder distance and determined the corresponding value of a. Since the angular 
diameter distance depends on a linearly for < z < 5, we define a as follows for the 
calculated distance D\ : 



a = (D A - D DR )/(D F - D DR ), 



(3-6) 



where Ddr is the limiting Dyer-Roeder distance with a = and Dp is the calculated 
Friedmann distance in the homogeneous case. This Dp is equal to the Dyer-Roeder 
distance with a = 1, except for small errors in numerical integrations. 

As a result of statistical analysis for this ray-shooting, we derived the average 
dumpiness parameter a, the dispersion a a , and the distribution (N(a)) of a. Here 
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Fig. 7. The percentage (WON(a)/N) of the 
distribution of a in bins with the inter- 
val Act = 0.4, for the parameter set A(20, 
0.005) in model S with (J? ,A ) = (1.0,0). 
Results for z — 0.5, 1, 2, 3, 4, 5 are denoted 
by dot-long dashed, dot-short dashed, long 
dashed, short dashed, dotted, solid lines, 
respectively 



Fig. 8. The percentage (100N(a)/N) of the 
distribution of a in bins with the inter- 
val Aa — 0.4, for the parameter set B(40, 
0.005) in model S with (1.0, 0). Lines have 
the same meaning as in Fig. 7. 



in order to study the frequency of ray pairs with a, we consider many bins with the 
interval Aa = 0.4 and the centers aj = 1.0 ± OAi (i = 0, 1,2, ...). The number of 
ray pairs with ai — Aa/2 < a < on + Aa/2 is N(cii) and the total number of ray 
pairs is N(= 2~2i^( a i))- I n Tables I ~ IV we show (a,a a ) for five parameters sets 
in model S, 01, 02 and L, respectively. In Figs. 7 ~ 18, we show the percentages of 
the distribution of a, that is, 100iV(a)/A^ for several parameter sets. The following 
types of statistical behavior are found from these tables and figures : 

(1) a is nearly equal to 1 in all models. However the individual values differ con- 
siderably from 1 for 9 < 1 arcsec, because the dispersion a a is large enough (for 
example a a ~ 0.5 and 1 for z = 1 — 2 in model S and low-density models, respec- 
tively). Of the four models the dispersions are smallest and largest in models S and 
01, respectively. Those in 02 and L are comparable. This situation is reflected also 
in the behaviors of the distribution of a, that is, around the point a = 1 it is most 
symmetric and asymmetric in models S and 01, respectively. 

(2) The dependence of a a on the radius of particles is large, a a decreases with an 
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Fig. 9. The percentage (WON(a)/N) of the 
distribution of a in bins with the inter- 
val Aa = 0.4, for the parameter set A(20, 
0.005) in model 01 with (0.2,0). Lines 
have the same meaning as in Fig. 7. 



Fig. 10. The percentage (W0N(a)/N) of the 
distribution of a in bins with the inter- 
val Aa — 0.4, for the parameter set 33(40, 
0.005) in model 01 with (0.2, 0). Lines 
have the same meaning as in Fig. 7. 



increase of r s , and when r s « 20/i kpc, the dispersion a a is too large to determine 
a. The symmetry of the distribution of a increases in the order of C, A, and B 
(for the same 9) in all models. This is because with an increase of r s the volume 
of empty regions included in the separation angles of ray pairs decreases and the 
inhomogeneous matter distribution in them is smoothed. 

(3) With the increase of the redshift z, the dispersion a a decreases and the symmetry 
of the distribution increases in all models. This is because with an increase of 
z the number of particles included within the separation angles increases and the 
gravitational effect on ray pairs is homogenized. 

(4) For 9 < 1 arcsec, the dependence of a a on the separation angle 9 is small, as 
can be found by comparing columns A and D of the tables. In low-density models 
the dependence is rather large for 9 >> 1 arcsec, as in column E. This is because 
in low-density models the quantity of matter included within the separation angle 
of ray pairs is much different in the cases of 9 < 1 and 9 » 1 arcsec, and the 
corresponding distance for 9 >> 1 arcsec is clearly close to the Friedmann distance, 
in contrast to the case of 9 < 1 arcsec. 

(5) In the case of a a > 0.8, the distribution of a is irregular and asymmetric around 
a = 1 and generally N(a) for a < 1 is larger than that for a > 1 (cf. Figs. 9 ~ 11, 
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Fig. 11. The percentage (WON(a)/N) of the 
distribution of a in bins with the inter- 
val Aa = 0.4, for the parameter set C(10, 
0.005) in model 01 with (0.2,0). Lines 
have the same meaning as in Fig. 7. 



Fig. 12. The percentage (W0N(a)/N) of the 
distribution of a in bins with the interval 
Aa = 0.4, for the parameter set E(20, 20) 
in model Ol with (0.2, 0). Lines have the 
same meaning as in Fig. 7. 



13 ~ 17). This behavior suggests that the distribution function of a is not Gaussian 
in these cases. For ray pairs with a a > 0.8, the z dependence of a also is irregular 
and may decrease slightly from 1 with an increase of z. For ray pairs with a a < 0.5, 
however, a clearly approaches 1 with an increase of z. 

(6) A positive cosmological constant has the role of decreasing a a and the asymmetry 
of the distribution of a (cf. Figs. 9 ~ 12 and Figs. 15 ~ 18). 

Figures 3 ~ 6 in the previous section show for z ~ 0.5 that the change in the 
Dyer-Roeder distance for the interval a = ~ 1 is very small. Accordingly, we 
can approximately use Da with a = 1 for z ~ 0.5, even if a a is larger than 1. For 
z > 1 and 6 < 1 arcsec, however, the change in Da for this interval of a is not so 
small that the large dispersion (a a ) brings a considerable fluctuation in the value of 
.Da- This means that for z > 1 and 9 < 1 arcsec, the realistic disntance may deviate 
considerably from the Friedmann distance in low-density models. 



§4. Concluding remarks 
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Fig. 13. The percentage (WON(a)/N) of the 
distribution of a in bins with the inter- 
val Aa = 0.4, for the parameter set A(20, 
0.005) in model 02 with (0.4,0). Lines 
have the same meaning as in Fig. 7. 



Fig. 14. The percentage (W0N(a)/N) of the 
distribution of a in bins with the inter- 
val Aa — 0.4, for the parameter set 33(40, 
0.005) in model 02 with (0.4, 0). Lines 
have the same meaning as in Fig. 7. 



Table I. The average dumpiness parameter a and its dispersion a a in model S with (!2o,Ao) = 
(1.0, 0). Columns A, B, C, D and E correspond to the parameter sets (20, 0.005), (40, 0.005), 
(10, 0.005), (20, 0.1) and (20, 20), respectively. 





A 




B 




C 




D 




E 




z 


Q 


a a 


a 


a a 


a 


a a 


a 


a a 


a 


a a 


0.5 


1.07 


0.93 


1.12 


0.78 


1.01 


1.31 


1.08 


0.84 


1.06 


0.76 


1 


1.01 


0.57 


1.05 


0.43 


0.98 


0.81 


1.02 


0.50 


1.01 


0.40 


2 


1.01 


0.36 


1.02 


0.28 


0.99 


0.54 


1.01 


0.33 


1.00 


0.26 


3 


1.00 


0.27 


1.02 


0.21 


0.99 


0.42 


1.01 


0.25 


1.00 


0.19 


4 


1.00 


0.23 


1.01 


0.18 


0.99 


0.35 


1.01 


0.21 


1.00 


0.16 


5 


1.00 


0.20 


1.01 


0.16 


0.99 


0.30 


1.00 


0.18 


1.00 


0.14 



By numerical ray-shooting in the iV-body-simulating clumpy cosmological mod- 
els, we studied the statistical behavior of the angular diameter distance D\ and 
determined the dumpiness parameter a by comparing it with the Friedmann dis- 
tance (a = 1) and the Dyer-Roeder distance (a > 0). The results show that the 
average value of a is nearly equal to 1 and the dispersion (a a ) decreases with an 
increase of z, but that for 9 < 1 arcsec, a a is ~ 0.5 in model S and ~ 1 in low-density 
models for z = 1 — 2. Hence the influence of the dumpiness is not small, because the 
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Fig. 15. The percentage (W0N(a)/N) of the 
distribution of a in bins with the inter- 
val Act = 0.4, for the parameter set A(20, 
0.005) in model L with (0.2,0.8). Lines 
have the same meaning as in Fig. 7. 



Fig. 16. The percentage (W0N(a)/N) of the 
distribution of a in bins with the inter- 
val Aa — 0.4, for the parameter set 33(40, 
0.005) in model L with (0.2, 0.8). Lines 
have the same meaning as in Fig. 7. 



The i 


iverage 


dumpiness parameter 


a anc 


its d 


ispersion 


(7a in 


model 


01 wit! 




A 




B 




C 




D 




E 




z 


a 




a 


a a 


a 


a a 


a 


a a 


a 


a a 


0.5 


1.01 


2.42 


1.23 


1.28 


1.07 


4.08 


0.99 


2.39 


0.91 


1.48 


1 


0.94 


1.43 


1.05 


0.86 


0.91 


2.38 


0.93 


1.42 


0.96 


1.08 


2 


1.00 


1.03 


1.02 


0.73 


0.95 


1.63 


0.99 


0.99 


1.00 


0.78 


3 


0.98 


0.89 


1.01 


0.63 


0.95 


1.37 


0.99 


0.87 


1.00 


0.66 


4 


0.98 


0.85 


0.99 


0.59 


0.93 


1.31 


0.98 


0.83 


1.00 


0.60 


5 


0.96 


0.80 


0.97 


0.55 


0.89 


1.26 


0.97 


0.78 


1.00 


0.54 



difference between distances with a = and 1 is not small for z > 1. For 9 » 1 arc- 
sec, a a is so small that the distance can be regarded approximately as the Friedmann 
distance. Therefore we can conclude for 6 < 1 arcsec that for rough or qualitative 
estimates of the lensing effect we can use the Friedmann angular diameter distance 
(a = 1), but for the quantitative analysis of cosmological gravitational lensing, some 
errors may result from using only it at the high-redshift stage (z > 1). 

For the behavior of the distribution (N(a)), we found that when a a < 0.5, N(a) 
is symmetric around a = 1, and when a a > 0.8, N(a) is irregular, and the symmetry 
is largest and smallest for S and 01 models, respectively. 
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Fig. 17. The percentage (WON(a)/N) of the 
distribution of a in bins with the inter- 
val Act = 0.4, for the parameter set C(10, 
0.005) in model L with (0.2,0.8). Lines 
have the same meaning as in Fig. 7. 



Fig. 18. The percentage (W0N(a)/N) of the 
distribution of a in bins with the interval 
Aa = 0.4, for the parameter set E(20, 20) 
in model L with (0.2, 0.8). Lines have the 
same meaning as in Fig. 7. 



Table III. The average dumpiness parameter a and its dispersion a a in model 02 with (0.4, 0). 

A B C D E 

z a a a a a a a a a a a a a a a 

~05 L00 L01 L00 079 099 L71 L00 097 L00 0.86 

1 1.02 0.89 1.00 0.65 1.03 1.50 1.01 0.86 1.00 0.65 

2 0.99 0.62 0.98 0.46 1.01 1.01 0.99 0.60 0.98 0.43 

3 0.99 0.51 0.99 0.37 0.99 0.80 1.00 0.50 0.99 0.34 

4 0.99 0.45 0.99 0.32 0.97 0.68 0.99 0.44 0.99 0.30 

5 0.99 0.40 0.99 0.29 0.97 0.63 0.99 0.38 0.99 0.27 



In the present cosmological models all particles were assumed to be lens objects 
with equal masses and radii. As the number density of particles is larger than 
that of visible standard galaxies, most particles are regarded as invisible galaxies 
or dark matter balls. It is a crucial problem in cosmological lensing to specify how 
strong lenses these invisible objects are. If the particles corresponding to only visible 
galaxies are regarded as lens objects, a a may be smaller than that in the above case 
A. If all visible galaxies are lens objects with r s = 20h~ 1 kpc and the other particles 
are weaker lenses with r s ~ 40/i _1 kpc, the resultant values for a and a a are between 
those in cases A and B. 
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Table IV. The average dumpiness parameter a and its dispersion a a in model L with (0.2, 0.8). 





A 




B 




C 




D 




E 




z 


a 


<J a 


a 


<J a 


a 


a a 


a 




Ol 


a a 


0.5 


1.01 


1.44 


1.03 


1.03 


1.00 


2.49 


1.02 


1.35 


0.97 


1.08 


1 


1.01 


1.01 


1.01 


0.64 


0.95 


1.78 


1.02 


0.96 


0.99 


0.60 


2 


1.00 


0.71 


1.00 


0.42 


0.94 


1.23 


1.00 


0.68 


0.99 


0.36 


3 


1.00 


0.53 


1.00 


0.31 


0.96 


0.92 


1.00 


0.51 


0.99 


0.26 


4 


1.00 


0.42 


1.00 


0.25 


0.95 


0.73 


1.00 


0.41 


0.99 


0.21 


5 


1.00 


0.36 


1.00 


0.21 


0.95 


0.63 


1.00 


0.35 


0.99 


0.17 



In the above averaging process, all light rays were taken into account. If we 
consider only weakly deflected light rays as contributing to weak lensing, the disper- 
sion a a will be a slightly smaller than the values in the above tables. However, the 
contribution of strong lensing to a a is small because of its small frequency. 

In this paper we treated the case in which m is on the order of the standard 
galaxy. If m is on the order of the rich-cluster mass, the influence of lensing on the 
distance is of course much larger and gives much larger dispersions. However, this 
situation is unrealistic. 
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